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We report on first experimental signatures for chaos-assisted tunneling in a two-dimensional 
annular billiard. Measurements of microwave spectra from a superconducting cavity with high fre- 
quency resolution are combined with electromagnetic field distributions experimentally determined 
from a normal conducting twin cavity with high spatial resolution to resolve eigenmodes with prop- 
erly identified quantum numbers. Distributions of so-called quasi-doublet splittings serve as basic 
observables for the tunneling between whispering gallery type modes localized to congruent, but 
distinct tori which are coupled weakly to irregular eigenstates associated with the chaotic region in 
phase space. 
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For two decades a new kind of tunneling mechanism 
produces great interest, since it demonstrates how the dy- 
namical features of a classical Hamiltonian system effect 
the behavior of its quantum counterpart . This so- 
called "dynamical tunneling" occurs whenever a discrete 
symmetry of the system leads to distinct but symmetry 
related parts of the underlying classical phase space. In 
contrast to the well-known barrier tunneling, dynamical 
tunneling only depends upon the probability for such a 
quantum particle, although classically forbidden, to leave 
certain regions of phase space and travel into others. This 
basically involves the coupling strength between distinct 
phase space regions. In the special case of two symme- 
try related regular regions separated by a chaotic area 
in a mixed phase space, semiclassical quantization yields 
pairs of quantum states which are localized to the cor- 
responding sets of congruent, but distinct tori. These 
so-called quasi-doublets show a very sensitive splitting 
behavior which depends upon the coupling to irregular 
eigenstates associated with the intermediate chaotic sea. 
This non-direct, enhanced coupling of regular eigenmodes 
via chaotic ones is what defines chaos-assisted tunneling 
iQ-^ in its original sense . 

The aim here is to demonstrate for the first time that 
chaos-assisted tunneling can be observed experimentally, 
even for a case where the size of the splitting is several 
orders of magnitude below the typical mean level spacing 
of the system. For this purpose we performed measure- 
ments on superconducting as well as normal conduct- 
ing microwave cavities constituting a special family of 
Bohigas' annular billiard |^-[lot]. This system has been 
proven in very extensive and certainly also very accurate 
computer simulations, especially in Refs. to be a 

paradigm for chaos-assisted tunneling and provides ac- 
cess for experimental investigation. 

The two-dimensional geometry of the annular billiard 
is defined by two circles of radius r resp. R, the latter 
being set to unity, and the center displacement or eccen- 



tricity (5, see left part of Fig.l. In the following, only 
the special one-parameter family r + 6 = 0.75 will be 
considered, since it provides all the features which are 
relevant for chaos-assisted tunneling: From the classical 
point of view the system shows a transition from inte- 
grable {6 = 0) to mixed behavior {d > 0), thus develop- 
ing a growing chaoticity with increasing S. Furthermore, 
the discrete reflection symmetry leads to congruent but 
classically distinct regions in phase space. 
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FIG. 1. The annular biUiard for S = 0.20 and r = 0.55 
(left part) together with the corresponding Poincare surface 
of section (right part). Beside a large chaotic sea with stable 
islands in the center the phase space clearly displays two sym- 
metry related but dynamically distinct regular coastal regions 
for l^l > 0.75. Two examples of horizontal lines correspond- 
ing to whispering gallery trajectories (clockwise as well as 
counterclockwise) are shown. 

To demonstrate this, Fig.l (right part) shows a typ- 
ical Poincare surface of section for the configuration 
{d = 0.20/r — 0.55). Here the area preserving Birkhoff 
coordinates L (point of impact on the the outer circle) 
and S — sin a (angular momentum of the billiard par- 
ticle) have been used. Beside a large chaotic sea with 
chains of stable islands in the center, both of which are 
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influenced by a change of the eccentricity 5, two symme- 
try related but distinct neutraUy stable coastal regions 
for \S\ > 0.75 can be observed. Per construction of the 
family r + 6 = 0.75 these regular regions are invariant 
under variations of 5, since the corresponding trajecto- 
ries do not hit the inner circle. As a consequence S is 
conserved, indicated by horizontal lines in the surface of 
section. Those lines correspond to two so-called whisper- 
ing gallery trajectories From this, the only difference 
between the two distinct regular regions is the sign of 5, 
i.e. the sense of motion for the propagating particle. 

The fundamental question now accounts for the quan- 
tum counterpart of the classically forbidden transport 
between the distinct coastal regions: dynamical tunnel- 
ing. Since the coupling between both regions crucially 
depends upon the topology and the size of the chaotic 
sea, the system is in particular adequate to study chaos- 
assisted tunneling. But what is the basic observable 
for the tunneling strength in the corresponding quan- 
tum system? To answer this it is very instructive to 
start with the integrable case [8 = 0). Solving the 
Schrodinger equation with Dirichlet boundary conditions 
leads to eigenvalues kn,m and eigenstates Vl/n.m, with the 
angular momentum quantum number n and the radial 
quantum number m. Due to EBK quantization the 
property S = n/kn.m is the quantum angular momentum 
which has to be compared with the classical S = sin a 
in order to find the location of a certain quantum state 
in the phase space. While in the classical system the 
reflection symmetry of the billiard leads to two distinct 
but related regular regions with opposed sense of motion 
for the propagating particle (i.e. the whispering gallery 
trajectories clockwise and counterclockwise), the corre- 
sponding quantum eigenstates are organized in doublets 
ioi 5 — Q with two parities, even and odd, respectively. 
However, continuously increasing the eccentricity 5 sys- 
tematically destroys this doublet structure, yielding sin- 
gulets for states with S = n/k right within the chaotic 
sea {S < 0.75) and quasi-doublets on the remaining reg- 
ular coast {S > 0.75). As in the case of the well-known 
double-well potential the very small splitting of 

those quasi-doublets is directly determined by the classi- 
cally forbidden tunneling, thus presenting a very effective 
observable for the hardly accessible tunneling strength. 
Since the location of a certain quasi-doublet on the reg- 
ular coast (defined by S* = n/k) as well as the transport 
features of the chaotic sea (defined by the eccentricity 6) 
have a direct impact on the splitting, its systematic inves- 
tigation allows the experimental study of chaos-assisted 
tunneling in the annular billiard. 

As in earlier studies (for an overview, see ||ll|]), we 
simulated the quantum billiard by means of a two- 
dimensional electromagnetic microwave resonator of the 
same shape (see l.h.s. of Fig. 1). The measurements were 
divided into two parts: Taking in total three different 
conflgurations of the family r + S — 0.75 (i.e. S = 0.10, 
0.15 and 0.20) we performed on one hand experiments 
with a superconducting Niobium resonator (scaled to 



i? = i m) at 4.2 K in order to measure quasi-doublet 
splittings within the frequency range up to / = 20 GHz. 
The very high quality factor of up to Q « lO*' allows a 
resolution of F// « 10~^, where F is equal to the full 
width at half maximum of a resonance. For demonstra- 
tion. Fig. 2 shows transmission spectra of the supercon- 
ducting resonator in the vicinity of 9 GHz. Beside sev- 
eral singulets exactly one quasi-doublet can be observed 
which shows a small but systematic displacement with 
the eccentricity S. In all cases the very small splitting of 
the quasi-doublet is clearly detectable. This is only due 
to the high frequency resolution of the superconducting 
resonator. In this context, however, it is important to 
note, that the position of the exciting antennas has to be 
choosen very carefully in order to minimize the perturba- 
tion on the whispering gallery type modes in the coastal 
region (see Fig.l) and thus not to influence the size of 
their physical quasi-doublet splitting. Using antennas 
right within the whispering gallery region of the billiard 
(see sketch on top of Fig. 2) always produce "false" split- 
tings even for the concentric system {6 = 0) with twofold 
degenerate states. We therefore used antennas in the 
shadow region of the inner circle also preserving the sym- 
metry of the whole geometry, cf. Fig. 2. 
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FIG. 2. Transmission spectra around 9 GHz with vary- 
ing eccentricity 5. Among several singulets exactly one 
quasi-doublet, slightly moving with 5, can be observed. In 
the zooming circles the abszissa is stretched by a factor of 50 
in order to visualize the quasi-doublet splitting. Spectra at 
eccentricites 5 < 0.10 could not be realized with the present 
set up of antenna locations. 

For a proper identification of the quantum numbers 
{n\m) of the modes associated with the quasi-doublets on 
the other hand we used a normal conducting Copper twin 
of the Niobium billiard cavity. There we measured the 
corresponding wavefunction resp. electromagnetic field 
distributions from which the relevant quantum numbers 
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n and m could be deduced, even if they are far from being 
"good quantum numbers" in the given eccentric systems 
which are non-integrable. This second part of the experi- 
ment was based on a field perturbation method originally 
introduced in accelerat or phys ics and used successfully in 
billiard research before ||l3| , |l4[] . According to Slater's the- 
orem a small metallic body inside the cavity locally 
interacts with the electromagnetic field in such a way that 
a frequency shift of the excited mode results from the 
compensation of the non-equilibrium between the totally 
stored electric and magnetic field energy. This frequency 
shift 

5/ = /o - / = /o (^aEo^ - bHo^ (1) 

with respect to the unperturbed mode (index 0) directly 
depends upon the superposition of the squared electric 
and magnetic fields, Eq and Hq, respectively. Since the 
quantum wavefunction is related to the electric field only, 
the magnetic component has to be removed by a proper 
choice of the geometry constants a and h in Eq.(l) by 
choosing needle-like bodies (1.84 mm in length and 1.00 
mm in diameter). Moving the body across the whole 
two-dimensional surface of the billiard with a spatial res- 
olution of about one tenth of a wavelength by means of a 
guiding magnet and detecting the frequency shift df at 
each position, finally provides the complete field distri- 
bution. Examples of those for the configuration 5 = 0.20 
in the vicinity of 9 GHz are plotted in the upper part of 
Fig.3. 

Of the three distributions only the middle one with the 
quantum numbers n = 18 (36 field maxima in the polar 
direction) and m = 1 (one field maximum in the radial 
direction) is characteristic for the modes which are local- 
ized in the whispering gallery region (see Fig. 2). Con- 
trary to this the distributions on the l.h.s. and on the 
r.h.s. show a totally different pattern. The parity of the 
distributions is determined in the following way: If there 
is maximum field strength on the line which defines the 
reflection symmetry of the billiard, positive parity can be 
assigned to the mode, likewise negative parity for zero 
field strength on this line. 

Underneath the squared electric field strength distri- 
butions in Fig.3, the corresponding transmission spec- 
trum taken at 300 K with the normal conducting Copper 
cavity is shown. The three broad resonances associated 
with the field distributions are much better resolved in 
the measurement at 4.2 K of the superconducting Nio- 
bium twin cavity. The small displacement in frequency 
of the resonances in the two measurements is due to me- 
chanical imperfections of each individual cavity and po- 
sitioning errors of the respective inner circles within the 
resonators. Mechanical uncertainties of order ±100 /um 
relative to the radius of the outer circle R =125 mm are 
sufficient to account for the observed displacements. 

The spectrum at 4.2 K in Fig.3 shows that the reso- 
nance magnified in the insert is in fact one of the expected 



quasi-doublets characteristic for chaos-assisted tunnel- 
ing. Naturally, this quasi-doublet is not resolved in the 
spectrum taken at room temperature, and the field dis- 
tribution in the upper part of Fig.3 proves that of the 
two modes corresponding to the doublet in the particu- 
lar case considered the one with negative parity is excited 
stronger than the other. 

Calculating the quantum angular momentum S = n/k 
(with k = 2nR/\ — 2TTRf/co, where f denotes the cen- 
troid frequency of the quasi-doublet and cq the speed of 
light) finally yields the position of the whispering gallery 
type mode on the corresponding classical surface of sec- 
tion, see r.h.s. of Fig.l. In the case of mode (18|1) one 
obtains S w 0.77 characteristic for a mode in the so-called 
"beach region" |^ defined by the borderline S — 0.75 be- 
tween the chaotic sea and the regular coast, respectively. 
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FIG. 3. Matching of field strength plots taken at 300 K 
with high spatial resolution and microwave spectrum taken at 
4.2 K with high frequency resolution. The combined normal 
conducting/superconducting setups allow to measure highly 
resolved quasi-doublets including quantum numbers (n|m) 
and parity. See text. 

This comparison demonstrates that the measurements 
combine the high spatial resolution of about A/10 for the 
normal conducting billiard with the high frequency res- 
olution of about 1/Q for the superconducting one, thus 
allowing a very effective classification of regular quasi- 
doublets as well as chaotic singulets in the range up to ap- 
proximately 14 GHz, where the splittings become smaller 
than the resonance widths of the superconducting res- 
onator. The difference in frequency between the peaks of 
each quasi-doublet were estimated through a non-linear 
fitting to "skew Lorentzians" (see Eq. (4) in |p^ ). 

In what follows, we only consider the family with 
quantum numbers (n|l), since it consists of some 30 re- 
solved and undoubtedly identified quasi-doublets within 
the measured frequency range. To uncover effects due to 
chaos-assisted tunneling, the splitting of a certain quasi- 
doublet has to be analyzed as a function of its corre- 
sponding position in the classical phase space. As men- 
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tioned above, this position might be expressed in terms 
of S* = n/k as the quantized analog of the classical angu- 
lar momentum S = sin a, directly representing the very 
location in the underlying surface of section. The result- 
ing curve, again for the configuration i5 = 0.20, is shown 
in the lower part of Fig. 4. 
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FIG. 4. Distribution of normalized splittings with respect 
to the position of a certain quasi-doublet in classical phase 
space. The error on the splitting is generally smaller than the 
size of the squares except for the smallest splitting observed 
where an error bar is placed on the point. Beside a smooth 
transition from states with a large splitting right within the 
chaotic sea {S < 0.75) towards states with a very small split- 
ting on the regular coast (S > 0.75) a local maximum occurs 
in the direct vicinity of the beach, representing a very impres- 
sive signature for chaos-assisted tunneling. 



Here, the distribution of normalized splittings |A///| 
shows a very smooth transition from chaotic states de- 
fined by large splittings right within the chaotic sea 
{S < 0.75) to regular quasi-doublets with very small 
splittings in the classical coastal region {S > 0.75). Also 
the measured field distributions show an increasing reg- 
ularity with growing S, as can be seen from the exam- 
ples in the upper part of Fig. 4. Especially mode (29|1) 
is hardly distinguishable from the corresponding concen- 
tric mode (with splitting zero) although the system is 
strongly eccentric. 

Besides this global behavior a first strong signature 
of chaos-assisted tunneling can be observed in the par- 
ticular shape of the splitting curve: In the direct vicin- 
ity of the beach at 5* = 0.75 the quasi-doublets show a 
locally enhanced splitting amplitude, thus indicating a 
very effective coupling between the regular coast and the 
chaotic sea. As described above, this corresponds to a 
locally enhanced tunneling strength in the beach region 
as theoretically predicted in M. 
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FIG. 5. Distribution of normalized splittings in the 
zoomed beach region for different eccentricities S. The local 
maximum is rebuilt by all measured quasi-doublets forming a 
systematically rising and a randomly falling part below and 
above 5* — 0.75, respectively. 

To evaluate the influence of the chaoticity of the sys- 
tem. Fig. 5 shows a direct comparison of the splittings for 
all eccentricities 6 in the vicinity of the beach at S* = 0.75. 
Note that the splittings for different S not only enhance 
the visibility of the maximum, they furthermore reveal 
an additional feature of chaos- assisted tunneling: While 
splittings on the rising left part of the maximum, i.e. 
within the chaotic sea {S < 0.75), are distributed quite 
systematically — e.g. the data points of (5 = 0.20 al- 
ways correspond to the largest splittings — shows the 
falling right part large fluctuations for a given eccentric- 
ity S. This effect is also theoretically predicted 1^,^,^ and 
accounts for the high rate of anti-crossings with chaotic 
modes for high angular momenta S. Thus on the r.h.s. 
of 5 = 0.75 the tunneling strength shows a very random 
dependence on the eccentricity S leading to strong fluctu- 
ations in the distribution of splittings. Finally, for even 
larger values of S the splitting amplitudes are of the or- 
der of the inverse quality factor. A/// sa 1/Q ss 10~^ 
defining the resolution limit of the present setup. 

In summary, we have presented first experimental sig- 
natures for chaos-assisted tunneling in a billiard. As the 
basic observable we have investigated the splittings of 
quasi-doublets with respect to their position in the clas- 
sical phase space and their dependence on the eccentric- 
ity d. A local maximum in the vicinity of the beach 
region with a systematically rising and randomly falling 
part has been found which directly reflects the enhanced 
tunneling strength at this critical location between the 
regular coast and the chaotic sea. In this context, a com- 
bined experimental setup using normal as well as super- 
conducting billiards has offered a very effective tool for 
measuring highly resolved quasi-doublets with properly 
identified quantum numbers. 
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